AiA065045 


OJ 

o 

C-5 
UJ 

f zd 


RADAR  TRACK  DATA  CORRELATION  OR 
REACHABLE  SETS  REVISITED: 

THE  REACHABLE  STATE 


U.S.  RRmr 

miSSILE 

RESEARCH 


Richard  E.  Dickson 
Systems  Simulation  Directorate 
Technology  Laboratory 


nno 


DEYELOPmEnT 

commRnD 

ttl  : ; 


ilHEm 


FEB  28  1979 


isinnsily 

A 


Redstone  Arsenal,  Alabama  35809 


OMI  FORM  1000,  1 APR  77 


Approved  for  public  release; 
distribution  unlimited. 


79  02  27  005 


J 


OISIK>trnON  INtTflUCTIONt 


DESTROY  THIS  REPORT  WHEN  IT  IS  NO  LONGER  NEEDED.  DO  NOT 
RETURN  IT  TO  THE  ORIGINATOR. 


OttCLAlMCR 

THE  FINDINGS  IN  THIS  REPORT  ARE  NOT  TO  BE  CONSTRUED  AS  AN 
OFFICIAL  DEPARTMENT  OF  THE  ARMY  POSITION  UNLESS  SO 
DESIGNATED  BY  OTHER  AUTHORIZED  DOCUMENTS. 


TRADE  NAMES 

USE  OF  TRADE  NAMES  OR  MANUFACTURERS  IN  THIS  REPORT  DOES 
NOT  CONSTITU  TE  AN  OFFICIAL  INDORSEMENT  OR  APPROVAL  OF  THE 
USE  OF  SUCH  COMMERCIAL  HARDWARE  OR  SOFTWARE. 


UNCLASS TVTRn 


SECURITY  classification  OF  THIS  PAGE  (Whan  Data  Bntarad) 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 

T-TS-eiv' 


2.  GOVT  ACCESSION  NO 


i4.  TITLE  (and  Subtltla) 


Radar  Track  Data  Ccrrelation  or  Reachable 
Sets  Revls'^ted:  The  Reachable  State,  j 


7.  AUTHORT*) 


Richard 


E. /Dickson  | 


•.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Commander 

US  Army  Missile  Research  and  Development  Command 
ATTN:  DRDMI-TD 

Redstone  Arsenal.  Alabama  


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Commander 


US  Army  Missile  Research  and  Development  Command 
ATTN:  DRDMI-TI 

Redstone  Arsenal ^ Alabama — a5Ma. 


U.  MONITORING  AGENCY  NAME  S AOORESSfi/  dlllatant  boat  CcntntUnt  Otllca) 


IS.  OISTRieUVION  STATEMENT  fa/ thi*  XapOfi; 

Cleared  for  Public  Release; 
Unlimited  Distribution 


READ  INSTRUCTIONS 
BEFORE  CC»(PLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  S PFPinEt-gnmrocn 


L 


Technical  ^epwt 


6.  PERFORMING  ORG.  REPORT  NUMBER 


T-78-61 


S.  CONTRACTOR  GRANT  NUMBERfaJ 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

13.  mumper  of  PAOtS~- 


28 


1!>.  SECURITY  CLASS,  (ol  thia  taport) 

UNCLASSIFIED 


15a.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 


nj  3J-  -I'm  I 


'7.  DISTRIBUTION  STATEMENT  (of  ttia  abattaef  antarad  fa  Slock  30,  ff  dfffarant  from  Raport) 


IS.  supplementary  '40TES 


IS. 


KE^  WORDS  (Contfaua  on  ravataa  afda  ff  nacaaaatr  and  fdantffr  hr  Mock  numbar) 


Aircraft  Tracking 
Radar  Track  Data  Correlation 
Reachable  Sets 
^Air  Traffic  Coi  trol 


/• 


d 


ABSTRACT  (CaoOaua  am  tararm  aS*  *f 


t fdantffr  brbtocli 


Given  an  aircraft's  position  and  velocity  at  time  and  the  limits  ou  the 
rate  at  which  the  aircraft  executes  horizontal  turns,  vhat  is  the  set  of  all 
possible  position  i*-  -Ime  t,  given  the  chan;e  in  heading  at  time  t and  tha£ 
the  speed  is  jncL;  , 

\ 


i the  chan;e  in  heading  at  time  t and  that 

79  02  27  005 


DD,!;sr7.i«3  EDITION  OF  t NOV  SS  IS  OBSOLETE 

3?5 


UNCLASSIFIED 


M. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Bkon  Data  Bntarad) 


I 


TABLE  OF  CONTENTS 

Page 


1.  Introduction 7 

2.  Reachable  Sets  Revisited 8 

3.  The  Reachable  State  14 

4.  Results  and  Conclusions 24 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1.  Extrapolated  flight  path 9 

2.  Dimensionless  plot  of  inner  and  outer  boundaries 

relative  to  the  linearly  extrapolated  position  11 

3.  Reachable  set  for  a standard  turn 12 

4.  Reachable  set  for  a standard  turn  and  a speed  of 

300  meters  per  second 13 

5.  Reachable  set  containing  the  reachable  subset  whose 

rate  of  turn  is  half  the  maximum  16 

6.  Reachable  set  for  theta  max  of2n 17 

7.  Reachable  set  for  theta  max  of  v 17 

8.  Reachable  set  for  theta  max  of  ^2  18 

9.  Reachable  set  for  theta  max  of  ^4  18 

10.  Reachable  state  and  substate  for  114°  turn  to  the  right 20 

11.  Continuous  turn  geometry 21 

12.  Initial  and  final  turn  geometry 21 

13.  Turn  geometries  combined 22 

14.  Reachable  state  for  114°  tium  to  the  right  and  a 

theta  max  of  2 n 24 


2 


SYMBOLS 


a Acceleration  normal  to  the  velocity  vector, 

A Initial  Position.  For  a given  change  in  heading,  theta 
C The  point  in  the  reachable  set  reached  by  a continuous  turn. 

I The  point  on  the  inner  boundary  reached  by  a turn  at  the  maximum  rate. 

O The  point  on  the  outer  boundary  reached  by  a turn  at  the  maximum  rate, 

r Radius  of  turn  for  the  maximum  rate  of  turn  (Equation  (5)). 

t 'fime. 

to  Initial  time. 

At,  T Time  increment  or  delay. 

Ti  Time- spent  in  turning,  o < Tj  <T. 

V Speed,  magnitude  of  the  veloci  ty. 

Position  relative  to  the  initial  position 

X Along 

Y Normal  to 

the  linearly  extrapolated  flight  path. 

Position  relative  to  the  linearly  extrapolated  position 

X Along 

y Normal  to 

the  linearly  extrapolated  flight  path. 
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For  a given  change  in  heading  the  slope  of  the  line  containing 
a All  points  on  all  porsible  inner  boundaries. 

)3  All  points  on  all  possible  outer  boundaries. 

y All  points  between  the  inner  outer  boundaries  of  the  reachable  set. 
0 Change  in  heading  or  course. 

0 Rate  of  change  in  heading,  rate  of  turn  (Equation  (6)). 

SUBSCRIPTS 

( )o  Outer  boundary. 

( )i  Inner  boundary. 

( )c  Continuous  turn. 

( )max  Maximum. 

( )min  Minimum. 
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1.  INTRODUCTSON 

In  the  1980*8,  tactical  air  (traffic) 
control  and  air  defense  systems  will 
exchange  air  track  data  via  digital 
commimication  links.  Where  radar 
coverages  overlap,  the  forwarding 
delays  in  the  communication  links 
may  cause  problems  in  correlating 
locally  .leld  air  track  data  with  that 
forwarded  from  a remote  reporting 
unit.  One  is  reminded  of  Segal’s  Law: 
“A  man  with  one  watch  knov'^t  what 
time  it  is.  A man  with  two  watches  is 
never  sure.  ” 

A point  which  should  be  made 
explicit  is  that  the  problem  is  one  of 
correlation  of  processed  data;  that  is, 
the  data  has  already  been  filtered.  If 
one  wished  to  filter  again,  it  would  be 
necessary  to  know  the  filters  which 
had  already  been  used  since  they 
would  become  a part  of  the  plant  or 
observation  matrix.  Such  an  approach 
does  not  appear  practical.  Since  the 
radar  track  data  are  estimates,  it 
should  be  possible  to  treat  them  in  a 
deterministic  fashion. 

If  the  only  “maneuver”  one  could 
expect  was  no  maneuver,  there  would 
be  “little”  difficulty  in  correlating 
tracks,  but,  unfortunately,  this  is  not 
the  case.  Though  not  the  only 
maneuver  one  would  expect,  no 
maneuver  is  probably  the  most  likely 
“maneuver”  since  most  of  an  aircraft’s 


flight  time  would  be  spent  in  level 
flight  at  constant  speed  and  heading. 
Of  course,  this  iriplies  the  aircraft  is 
unalerted  or  cooperative. 

The  next  most  likely  maneuver 
would  seem  to  be  a change  in  coiurse  or 
heading,  that  is,  a single  turn  in  the 
horizontal  plane. 

The  above  maneuvers  are  rather 
“simple”  and  are  those  one  would 
expect  of  commercial  aircraft  or 
tactical  aircraft  in  formation.  Tactical 
aircraft  would  additionally  execute 
maneuvers  such  as  the  split  “S”, 
Immelman  turns,  loops,  diving  spirals, 
breaks,  diving  turns  and  the  like  when 
conditions  warrant. 

Hoares’  Law  of  Large  Problems: 
“Inside  every  large  problem  is  a small 
problem  struggling  to  get  out.”  In  the 
following,  the  emphasis  will  be  on 
single  horizontal  turns  at  constant 
speed.  As  in  Reference  2,  the  approach 
will  be  geometric  “because  it  provides 
more  insight  than  other  alternatives,” 
and,  also,  at  times  heuristic. 

2.  REACHABLE  SETS 
REVISITED 

In  Reference  1 and  2 the  question, 
“Given  the  aircraft  position  and 
velocity  at  time  t — A t and  the  limits  on 
the  ability  of  the  aircraft  to  maneuver 
and/or  change  speed,  what  is  the  set  of 
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all  possible  aircraft  positions  at  time 
t?”  is  addressed.  The  set  of  all  possible 
positions  was  called  the  reachable  set 
(or  set  of  attainability)  and  was 
described  by  a set  of  equations  for  the 
inner  and  outer  boundaries.  Only  the 
assumptions  and  results  are  presented 
here  since  the  derivation  is  readily 
available  [1]. 

The  following  assumptions  were 
made  concerning  the  motion  of  the 
aircraft  [1]: 

• **Any  lateral  ma'zeuver  in  the  xy~ 
plane  is  due  to  a constant-rate 
turn  at  constant  speed. 

• Only  one  turn  maneuver  is 
allowed  during  the  time 
inc.  ement  A t. 

• The  turning  maneuver  can  begin 
or  end  at  any  time  during  the  time 
increment  At. 

• Speed  variations  and  turning 
maneuvers  are  assumed  to  be 
independent  of  each  other. 

• Any  speed  change  is  assumed  to 
take  place  at  the  beginning  of  the 
time  increment  a t,  and  the 
speed  is  assumed  to  be  constant 
throughout  this  time  increment.” 

The  outer  boundary  of  the  reachable 
set  is  determined  by  initially  turning  to 
a new  heading  and  then  flying  out 
along  that  heading.  This  jaximizes 
the  departure  from  the  initial  position 
along  a given  heading  relative  to  the 
initial  position  [1].  Of  course,  flight 


along  the  initial  heading  maximizes 
the  departure  from  the  initial  position. 
The  inner  boundary  for  a single  turn  is 
determined  by  flying  along  the  initial 
heading  till  the  last  possible  moment 
before  executing  a change  in  heading; 
that  is,  the  turn  is  at  the  end  rather 
than  the  b^pnning. 

The  outer  boundary  of  the  reachable 
set  is  given  by  [1]. 

r sin  0 + (vAt-r*6)  cos  6(1) 

Yq  = ± [r  (1-cos  0)  + 

(v At-r* 0) sin  0] 

and  the  inner  boundary  by  [1]. 

= r sin  0 + (vAt-r*0)  (3) 
= + r (1  - cos  0)  (4) 

Since 

(5) 

and 


The  above  equations,  (I),  (2),  (3)  and 
(4),  may  be  rewritten 

2 /3T.\ 

sin(^j+v(T-Tp  cos^-^j 

(7) 
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Yi  = ±[r  (^)] 


It  should  be  noted  that  the  inner  and 
outer  boundaries  coincide  at  two 
points,  when  Ti  = T and  Ti  = 0. 

Since  the  most  likely  maneuver  may 
be  “no  maneuver”,  it  would  seem 
reasonable  to  relocate  the  origin  at  the 
linearly  extrapolated  position  of  the 
aircraft. 

Relative  to  the  linearly  extrapolated 
position,  the  outer  and  inner 
boundaries  are 


2 /aT,\ 

= I”  ( V ) \ V / 

(11) 


and 


respectively.  The  inner  boundary  is  no 
longer  a function  of  the  total  delay,  T, 
but  only  of  the  time  spent  in  turning, 
T^,  that  is,  the  inner  boundary 
comcide  for  all  time  delays. 

Perhaps  a gtdanken  experiment  will 
help  clarify  this  point.  Consider  a 
string  of  length  vT  and  a circle  whose 
radius  is  the  radius  of  turn, 

2 

^ - r (15) 
Let  the  string  be  attached  to  the  circle 
and  the  point  of  attachment  be  the 
tangent  point  between  the  circle  and 
the  linearly  extrapolated  flight  path 
(see  Figure  !)■ 


The  other  end  of  the  string  is 
attached  to  the  initial  position  of  the 
aircraft,  that  is,  the  position  at  time  t^. 
As  the  circle  rolls  (without  slipping) 
back  along  the  extrapolated  flight 
path,  the  string  is  wound  up  (without 
stretching)  on  the  circle.  The  locus  of 
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the  point  of  the  string’s  attachment  to 
the  circle  describes  a cycloid  [3],  the 
inner  boundary  of  the  reachable  set. 
When  the  string  is  all  wound  up,  the 
circle  will  stop  rolling  and  will  be 
tangent  to  the  linearly  extrapolated 
flight  path  at  the  aircraft’s  initial 
position.  The  other  half  of  the  inner 
boundary  would  be  constructed  in  a 
similar  fashion  on  the  other  side  of  the 
linearly  extrapolated  flight  path. 


y 0 

o max 

V 


= ( 1 - cos  0) 

+ (0  - 0)  sin  0 

max 


(18) 


and 


X.  0 
1 max 


= sin  0-0 


(19) 


• (1  - cos  0)  (20) 


Now  disconnect  the  string  from  the 
circle,  holding  the  circle  fixed  and 
keeping  tension  on  the  string,  unwind 
the  string  from  the  circle.  The  involutef 
generated  by  the  end  of  the  string  is  the 
outer  boundary  of  the  reachable  set 
[2]. 


respectively,  where 

« 

0 =0  T , /21> 

max  max 

and 

0=0  T,  . (22) 

max  1 


If,  for  some  reason,  one  knew  the 
linearly  extrapolated  position  but  not 
the  time  delay,  one  could  construct  the 
inner  boundary,  but  the  outer  bound- 
ary and  where  it  met  the  inner 
boundary  woul  i not  be  known. 

Since 

0 = -2.  (16) 
V 

The  equations  for  the  inner  and  outer 
boundaries  may  be  written 

_^oJmax  ^ e 

V 

+ (0  - 0)  cos  0—0 

max  / '-uo 

(17) 


Figure  2 is  a dimensionless  plot  of 
the  reachable  set;  actually  the 
dimensions  are  radians.  It  shows  the 
common  inner  boundary  and  various 
outer  boundaries.  In  Figure  3 the  total 
time  delays  for  the  outer  boundaries 
were  computed  using  Equation  (21)  and 
a standard  rate  of  turn  which  is  three 

degrees  per  second.  The  axes  now  also 
have  the  dimensions  of  time.  In  Figure 
4 the  speed  was  chosen  to  be  300  meters 
per  second  (mach  0.91, 583  knots).  One 
multiplies  the  speed  by  the  times  on  the 
axes  to  obtain  the  displacements.  Of 
course,  one  may  choose  any  rate  of 
turn,  and  give  a rate  of  turn,  any  speed. 


twhose  evolute  is  the  circle  [3]. 
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max 


V 


Figure  2.  Dimensionless  plot  of  inner  and  outer  boundaries  relative  to  the  linearly 
extrapolated  position. 


^ (sec) 


Figure  3.  Reachable  set  for  a standard  turn. 


Reachable  set  for  a standanJ  turn  and  a speed  of  300  nt/s.  millii 


To  construct  the  reachable  set,  it  was 
necessary  to  know  the  velocity  (speed 
and  heading)  and  position  at  the 
initial  time,  to  • addition  to  the 

position,  the  velocity  is  known  at  time 
t,  how  might  this  add’*ional 
information  be  used?  Obviously  one 
would  compare  the  speeds  to  determine 
if  they  were  equal,  ana,  if  equal,  is  the 
change  in  beading,  if  any,  within  the 
capabib'ties  of  the  aircraft? 


0 = arctan  (y/x)  (27) 

which  is  not  unexpected. 

What  is  the  slope,  r i of  the  line  con- 
necting points  on  the  inner  and  outer 
boundaries  for  a single  turn  that  have 
the  same  change  in  heading? 

tan  Y = 

yo  - y±  (28) 


3.  THE  REACHABLE  STATE 


~ cos  9-1 
sin  0 


(29) 


Given  the  aircraft’s  position,  speed 
and  heading  at  time  t^  and  the  limits 

on  the  rate  at  which  the  aircraft 
executes  horizontal  turns  and  that  the 
speed  is  unchanged,  what  is  the  set  of 
all  possible  positions  at  time  t given 
the  heading  at  time  t? 

Before  the  turn 

^ (23) 

y = o (24) 

but  after  a single  turn 

= V cos  Q (25) 

- fl  (26) 

y = V sin  0 ' ' 

and,  it  follows  that. 


= - tan  (0/2)  (30) 

Therefore,  the  magnitude  of  the  slope 
7 is  half  the  change  in  heading. 

From  another  gedanken  experiment 
it  will  become  apparent  that,  for  a 
given  rate  of  turn,  the  change  in 
headingis  the  same  at  all  points  along 
the  slope. 

Consider  the  case  where  the  rate  of 
turn  is  the  maximum  rate.  As  the  circle 
which  generated  the  involute  slides 
(without  rolling)  up  the  extrapolated 
flight  path,  the  end  of  the  string  moves 
linearly  from  the  outer  boundary  to  che 
inner  boundary.  Linearly  since  the 
amount  of  string  on  the  (arc  of  the) 
circle  remains  constant  and  any  string 
added  between  the  fixed  end  (initial 


I 


! 

I 

t 
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position)  and  the  circle  along  the 
extrapolated  flight  path  must  be  lost 
between  the  circle  and  the  free  end. 

The  reachable  set  for  any  turning 
rate  less  than  the  maximum  is  a subset 
of  (is  contained  in)  the  reachable  set  of 
interest  {Figure  5).  Each  of  the  possible 
subsets  would  contain  a similar  line 

segment  of  slope  r between  its  inner 
and  outer  boundaries.  The  combina- 
tion of  all  such  line  segments  would 
define  an  area  within  which  any  point 
in  the  reachable  set  with  a change  in 
heading  of  theta,  0,  would  lie.  The  line 
s^pnent  for  the  maximum  turning  rate 
discussed  above  defines  one  boundary 
of  the  area  in  question. 

Consider  the  slope  of  points  on  the 
inner  boundaries  of  all  reachable 
subsetst  for  a given  change  in  head- 
ing. Dividing  Equation  (13)  by 
Equation  (14) 

(9*^)  - (62)  sin  0-9 

: r-=  (31) 

and  taking  the  arctangent 


a = 


arctan 


' sin  0-0 


1 - cos  0 


gives  the  desired  slope.  Equation  (32) 
does  not  appear  very  pleasant  but  a 
“small  angle”  approximation  makes  it 
less  formable. 


t for  all  possible  turning  rates  S Q max. 


a 


1 


2 

- 9_  + 
2! 


(33) 


(34) 


A check  on  the  accuracy  of  the 
approximation  for  various  changes  in 
heading  will  be  found  in  Table  1. 

The  “small  angle”  approximation  is 
quite  excellent,  that  is,  “good  enough 
for  government  work.”  Beyond  180° 
the  approximation  breaks  down 

seriously.  Tactical  computers  will  be 
pressed  for  CPU  time  and  memory  and 
approximations  can  reduce  the 
computational  burden. 

The  slope  7 is  the  slope  of  the  line 
segment  between  points  on  the  inner 
and  outer  boundary  that  have  the 
same  change  in  heading.  The  slope  a 
is  the  slope  of  the  line  segment  between 
the  inner  boundary  and  the  point 
reached  by  a continuous  turn  whose 
change  in  heading  is  theta.  The  mini- 
mum turning  rate  for  a change  in 
heading  theta  is 

®min  ~ f (35) 

It  follows  from  Equation  (13)  and 
(14)  that  the  point  reached  by 
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continuously  turning  at  the  minimum 
turning  rate  is 


x-  = — (sin  e - 0)  (36) 

9 


From  Equations  (11)  and  (12) 

-yo(®2>  ‘ 


yc  - — (1  - cos  6) 
0 


8ln0  -00080 

(1  - CO80)  - 08in0 


If  nondimensional  plots  are  desired,  and  the  slope  is 


^max  = ^max  (sin  0-0) 


Vc  0 0 

niax  = max  (1  - cos  0) 


Figures  6 through  9 are  nondimen- 
sional plots  like  Figure  1 with 
“continuous  turns”  added.  One  would 
follow  the  same  procedure  used  in 
Figures  3 and  4 to  inti  >duce  the  desired 
rate  of  turn  and  speed. 

Equations  (38)  and  (39)  v/ill  generate 
the  curve  containing  all  points  reached 
by  a continuous  turn  within  the  reach- 
able set. 

What  remains  to  be  determined  is 
the  slope  of  the  line  containing  all 
possible  points  on  the  outer  boundaries 
of  all  possible  reachable  subsets  for  a 
given  change  in  heading  and  time 
delay  for  all  possible  turning  rates 
(Figure  10). 


Q _ « |3ln9  - 9cos9  *1 

8 -arctan|^^^  -co80)  - 08ln0j  (41) 

If  Equation  (32)  is  not  very  pleasant, 
Equation  (41)  is  positively  unpleasant. 
“Small  angle”  approximation  yeilds 


Obviously 


0+^  = 0 
3 3 


and  one  might  suspect 

a + B = 0,  (44) 

which  is  indeed  true  as  will  be 
demonstrated. 

First  consider  the  case  of  a 
continuous  turn,  Figure  11. 

AC  = 2 r sin  (6/2)  (45) 

X = (2  r sin  (9/2))  cos  0/2  (46) 

y = (2  r sin  (0/2))  sin  9/2  (47) 


Also  points  on  the  inner  and  outer 
boundary  that  have  the  same  change 
in  heading  are  the  same  distance  from 
the  initial  position  (Figure  12). 

The  arc  B I is  equal  to  the  arc  A D 
and  A B is  equal  D 0,  therefore  A I is 
equal  A O. 

It  follows  that  D I equals  D O, 
therefore  I M is  equal  M C,  and  I C is 
equal  C O, ... , and,  finally, 

i.a-B-f  <«) 

2 2 

which  i mplifies  to 


maximum  turning  rate  of  three  degrees 
per  second 

T < 60s 

will  not  result  in  difficulties.  For  longer 
delays  all  possibilities  should  be 
checked. 

Since  sufficient  time  must  have 
elapsed  for  the  change  in  heading  to 
have  occurred,  one  should  first  check 
to  insure  that 

9 =9  T 

max  max 

is  not  exceeded.  This  test  will  eliminate 
points  which  are  obviously  not 
members  of  the  reachable  state. 


a + ti  = e • <49) 

Equation  (44)  would  be  useful  if  one 
had  computed  a using  Equation  (32) 
since  one  may  determine  0 using 
Equation  (44)  instead  of  Equation  (41). 

A word  of  caution;  a change  of 
heading,  theta,  could  result  from  a 
“simple”  turn  or  a turn  of 

±(2  n IT  + 9)  , n > o . 

if  time  allows.  For  example,  a change 
of  heading  of  ninety  degrees  to  the 
right  could  have  resulted  from  a ninety 
degree  turn  to  the  right,  a two  hundred 
and  seventy  degree  turn  to  the  left,  et 
cetera.  For 

9 T < TT. 
max 

there  is  no  difficulty,  which  is  also  the 
practical  limit  of  the  small  angle 
approximations  given  above.  For  a 


The  reachable  state  for  a one 
hundred  and  fourteen  degree  turn  to 
the  right  and  theta  max  of  two  pi  is 
illustrated  in  Figure  14.  For  clarity  the 
state  for  two  hundred  and  forty-six 
degrees  to  the  left,  360°  - 114°,  is  not 
shown.  How  the  reachable  state  may 
be  implemented  is  as  follows. 

If  the  point  P (X,  Y)  lies  in  the  reach- 
able state 


tan 


\y  - y. ' 


tan  7 (50) 


and 


> tan  3 


(51) 


7 . . 


23 


2tt  IT  IT  2ir 


Figure  14.  Reachable  state  foi  114**  turn  to  the  right  and  a theta  max  of  2 ir  . 


Substituting,  Equation  (50)  and  (51) 
become 


respectively. 


4.  RESULTS  AND 
CONCLUSIONS 

It  was  found  advantageous  to  locate 
the  origin  of  the  reachable  set  at  the 
linearly  extrapolated  position;  the 
inner  boundary  is  then  common  to  all 
sets  for  a single  turn.  This  then  led  to 
the  reachable  state  where  the  change 
in  heading  is  used  to  further  reduce  the 
allowable  area  in  which  the  position 
may  lie.  The  triangle  of  the  reachable 
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state  sweeps  out  the  reachable  set  as 
one  varies  the  change  in  heading  from 
zero  to  the  maximum  possible  for  the 
time  delay. 

A point  which  should  be  made 
explicit  is  that  the  reachable 
sets/states  not  only  answer  the 
question  as  to  what  the  aircraft  could 
attain  but  how  its  present  position 
could  have  been  attained,  that  is, 
where  it  could  have  come  from.  Since 
one  may  be  correlating  track  data 
which  are  members  of  sequences,  those 
closest  together  in  time,  whether  in  the 
relative  “past”  or  “future”,  should  be 
checked. 

A caveat,  the  only  conclusion  one 
may  draw  if  a point  fails  the  reachable 
state  criteria,  is  that  it  does  not 
correlate  with  an  aircraft  making  a 
single  horizontal  turn  at  constant 
speed. 

It  appears  that  multiple  turns  could 
be  addressed  in  a similar  manner.  The 
outer  boundary  would  remain  the 
same  since  it  is  a maximum  [1]. 
However,  for  two  opposing  turns  the 
inner  boundary  would  be  composed  of 
epicycloids  at  the  initial  position.  One 
would  roll  one  circle  of  minimum 
radius  around  another  on  the  other 
side  of  the  linearly  extrapolated  flight 
path  till  the  point  reached  the  linearly 
extrapolated  flight  path  and  then 
repeat  the  process  from  the  other  side. 


Heuristically  this  appears  to  be  the 
inner  boundary  for  multiple  turns  but 
a proof  has  not  been  developed  at  this 
time.  The  reachable  state  for  two 
opposing  turns  resembles  the 
reachable  set  for  a single  turn  with  the 
radius  of  turn  doubled  [4].  The  reach- 
able set  for  double/multiple  turns 
would  have  uses  in  determining 
acquisition  cells.  The  equations  are  : 

Xi0 

— — » 2sin6  - sin  (6max-TT-20) 

(54) 

Yi® 

— — = 1 - 2cos6  - cos(0max-iT-28) 

(55) 

The  effects  of  uncertainties  remain 
to  be  incorporated.  The  method 
suggested  in  Reference  1 would  be 
rather  pessimistic,  but  for 
determination  of  acquisition  cells  this 
simplicity  may  be  acceptable. 

The  gedanken  experiment  is  worth 
performing.  The  Kenner  Spirograph® 
is  readily  available  or  consult 
Reference  3 for  other  methods  — or 
design  your  own;  a piece  of  string  ... 

The  reachable  sets/states  for  single 
turns  are  a possible  means  of 
correlation  when  the  length  of  the  time 
delays  warrant.  At  the  very  least,  they 
graphically  illustrate  the  effects  of 
time  delays  on  track  data.  In  addition 
to  fixed  wing  aircraft  there  would  also 
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be  applications  to  tracking  surface 
ships.  Though  motivated  by  the  needs 
of  tactical  air  defense  and  control,  the 


potential  for  such  techniques  in 
commercial  air  traffic  control  is 
obvious. 


t 
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